Intrinsic Superconductivity at 25 K in Highly Oriented Pyrolytic Graphite 
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High resolution magnetoresistance data in highly oriented pyrolytic graphite thin samples man- 
ifest non-homogenous superconductivity with critical temperature T c ~ 25 K. These data exhibit: 
i) hysteretic loops of resistance versus magnetic field similar to Josephson-coupled grains, ii) quan- 
tum Andreev's resonances and iii) absence of the Schubnikov-de Haas oscillations. The results 
indicate that graphite is a system with non-percolative superconducting domains immersed in a 
semiconducting-like matrix. As possible origin of the superconductivity in graphite we discuss 
interior-gap superconductivity when two very different electronic masses are present. 

PACS numbers: 74.10.+v,74.45.+c,74.78.Na 



The standard way to ascribe superconductivity to ma- 
terials is by observing the screening of an external ap- 
plied magnetic field, the Meissner effect, below a critical 
field B c i and, although less important from the physi- 
cal point of view, by measuring the drop of resistance to 
practically zero below a critical temperature T c . These 
phenomena are observed for percolative or homogenous 
superconductors where a macroscopic wave function of 
the Cooper pairs exists It is well known that in in- 
homogeneous superconducting samples, as for example 
the well-known ceramic high T c oxides, sometimes su- 
perconductivity does not percolate, then the resistance 
does not drop to zero and the Meissner effect is small. 
In this case the criteria to assign non-percolative inho- 
mogeneous superconductivity to a material is much less 
obvious. In addition, we would like to discuss here a su- 
perconducting high-T c material with a very low density 
of free electrons or quasiparticles n < 10 18 cm~ 3 , with 
very different effective masses to*. We think that this is 
the case of highly oriented pyrolytic graphite (HOPG), 
the material studied in this work. 

Untreated HOPG samples manifest large electronic 
mean free path and Fermi wavelength of order of microns 
[2| • On the other hand the same samples reveal that the 
surface is not an equipotential with metallic and insu- 
lating regions that can move 0, 0] . It seems clear that 
the view of graphite as a more or less ordered, homoge- 
neous system and with a homogeneous density of carriers 
cannot be hold and it does not represent the interesting 
piece of the physics of HOPG. Although resistance R(T) 
data can be fitted, in some cases, with an homogeneous 
two band model (TBM) using two mobilities and two car- 
rier concentrations (all temperature dependent parame- 
ters) Q , there are other observations as a function of the 
applied magnetic field reported here that cannot be ex- 
plained within this model. In this work we treat HOPG 
as a non-uniform electronic system and as such it will be 



discussed. 

To aboard this hard problem we have obtained over 
10 6 high resolution magnetoresistance (MR) data points 
in a range of temperatures. These data exhibit: (i) ir- 
reversible hysteretic loops of resistance versus magnetic 
field similar to those observed in granular superconduc- 
tors with Josephson-coupled grains 0, 0] that can be 
assigned to superconducting fluxons, (ii) quantum An- 
dreev's resonances in the MR Q and (iii) absence of 
Schubnikov-de Haas (SdH) oscillations. The experimen- 
tal data indicate the existence of energy gaps at the Fermi 
level and that HOPG is a non-percolative superconductor 
with "granular" domains immersed in a semiconductor- 
like matrix. The origin of the superconductivity in 
graphite may be assigned to interior-gap superconduc- 
tivity that predicts a gapless stability when two different 
masses are present^a problem that has been discussed by 
Liu and Wilczeck [9( . 

The high-resolution, low-noise four-wires MR measure- 
ments have been performed by AC technique (Linear 
Research LR-700 Bridge with 8 channels LR-720 multi- 
plexer) with ppm resolution and in some cases also with 
a DC technique (Keithley 2182 with 2001 Nanovoltmeter 
and Keithley 6221 current source). The temperature sta- 
bility achieved was ~ 0.1 mK and the magnetic field, 
always applied normal to the graphene planes, was mea- 
sured by a Hall sensor just before and after measuring 
the resistance, and located at the same sample holder 
inside a superconducting-coil magnetocryostat. We used 
currents between 1 . . . 100 /iA. 

To start with our strategy we have prepared different 
samples of HPOG that just differ in its ordering and size 
and they exhibit apparently different behaviors with T. 
Figure Q] shows R(T) for the samples indicated in the fig- 
ure caption. Usually one tends to fit these curves with 
the TBM. In particular the R(T) of sample (3) can be fit- 
ted approximately. However, carriers in HOPG have two 
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different masses and one of them is practically zero corre- 
sponding to Dirac electrons [10]. Furthermore, there are 
other important aspects described below that undoubt- 
edly cannot be put into accord with the TBM. We con- 
centrate in the very thin and micrometer small sample 
because it should have less number of fluctuating domains 
and this should provide more clear superconducting- 
related effects. Note that this sample shows a semicon- 
ducting like behavior that levels off at T ~ 25 K; its in- 
plane resistivity p a b( 10 K) ~ (50 ± 10) /iilcm is similar 
to the one of sample (1) from which it has been obtained 
by careful exfoliation. 

Figure[2fa) shows the MR of sample (3) at 4 K in detail 
and in the region 4 T to 8 T with larger resolution using 
a magnetic field step of ~ 1 Oe. The first surprise is that 
the MR is very small compared with the MR of larger 
samples of HOPG. In these samples the ordinary MR of 
HOPG between T and 8 T is - 10000% while in the 
small sample measured here is only < 300%. This differ- 
ence is discussed in Ref. H. In addition, SdH oscillations 
are absent in sample (3) (in other samples of similar size 
we measured they appear very weak). This might imply 
that the Fermi level lies in a gap. Notice that we decided 
to perform experiments with very small field increment. 
This was not done accidentally. The reason is that we 
expected to have weak quantum oscillation resonances 
- compared with the classical SdH oscillations - due to 
the small number of potential fluctuations (note that the 
sample is small, of the order or smaller than the mean free 
path and Fermi wave length) and these fluctuations will 
induce an oscillating transmissivity through the potential 
wells. These quantum oscillations were proposed theoret- 
ically to interpret observed structures that were over seen 
or consider noise in graphene samples Q . And of course 
the sample of Fig. [2] shows the expected quantum oscil- 
lations. These quantum oscillations have a two period 
spectrum indicating that in the sample one has at least 
two characteristic potential wells. Figure [Hh) shows the 
oscillation amplitude of the two harmonics (see also the 
inset) as a function of T, which remain constant below 
10 K and vanish at a critical temperature T c ~ 25 K. 

We claim that these oscillations, given their small am- 
plitude of ~ 100 nV. . . 400 nV (much smaller than the 
corresponding values in temperature for the used range 
T ^ 2 K) are due to the interference of wave functions 
that suffer Andreev's reflections at the potential walls 
matching low-gap semiconducting with superconducting 
regions. From the period of the oscillations in field we 
can estimate that there are superconducting "granular" 
domains of size around 1 /im separated by small-gap 
semiconducting matrix of similar size, which couples the 
superconducting grains. If this picture is realized one 
expects to sec pinning and dissipation effects due to flux- 
ons, as discussed by Ji et al. in Ref. @, with circumvent 
superconducting currents between the superconducting 
grains through the semiconducting regions. One may 




0.0 



massive electrons 
Dirac Fermions ^ 
Wavevector k 



50 100 150 200 250 
Temperature T(K) 



FIG. 1: Normalized resistance as a function of temperature 
at zero applied field for three HOPG samples. Samples (1) 
and (3) were obtained from the same HOPG grade ZYA (0.4° 
rocking curve width) bulk sample from Advanced Ceramics. 
Particle induced x-ray emission indicates impurity levels of 
metallic elements below 5 [ig/g with exception of V (16 flg/g). 
Sample (1) (#(275 K) = 72 mf2) was ~ 10 /im thick and 
4.4 mm length. Sample (3) (#(275 K) = 16 ft) was 12 nm 
thick and ~ 30 /im length and width, with ~ 9 /tm distance 
between nearby electrodes, see inset in Fig. [2ja). The Pd- 
electrodes (to avoid Schottky barriers) were prepared using 
conventional electron lithography. The sample (2) (_R(275 K) 
= 6.5 mf2) was obtained from HOPG grade ZYC bulk sample 
(3.5°). The inset shows schematically the energy dispersion 
relations for two carriers in graphite, massive and massless 
(Dirac fermions). Following Ref. @, the instability region lies 
around the Fermi wavevector of the light particles. 



argue against the physical ground of the model we are 
proposing: how is it possible that superconducting pairs 
can be kept in a micron-size semiconducting-like regions 
connecting the superconducting ones? This should not be 
a problem. By using nano-fabricated constrictions and 
measuring the transition from ohmic to ballistic trans- 
port we have observed that the mean free path of the 
carriers in HOPG at 10 K is > 10 fim. Therefore, it 
should be perfectly possible that the pairs travel ~ 1 /im 
distance without breaking out. In other words the prox- 
imity effect in graphite may extend to microns. 
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FIG. 2: (a) Resistance of sample (3) between two adjacent 
voltage electrodes as a function of magnetic field. A close in- 
spection of the MR of this sample at fields above 0.5 T reveals 
an anomalous behavior, namely the MR oscillates. The os- 
cillations shown in the insets were obtained after subtracting 
a quadratic field dependence around 5.3 T and 6.4 T. These 
small-field-period oscillations in the resistance are superposed 
to oscillations of larger amplitude and field period, see inset 
in (b). Further measurements indicate that the overall shape, 
field positions and period are independent of the field sweep- 
ing rate, field step and field sweep direction. The oscillations 
are observed at low as well as high fields, as expected be- 
cause the slope of R vs. B does not depend appreciable with 
field [IJ. Different periods as well as oscillation amplitudes are 
observed for other samples. The inset shows an optical micro- 
scope picture of the sample with the Pd-electrodes. Count- 
ing clockwise from input current electrode 1 at the right, the 
data shown were taken between electrodes 3 and 4 (Ch.2). (b) 
Temperature dependence of the voltage amplitude of the two 
oscillations taken from the Fourier fit, see inset. The contin- 
uous lines are a guide. The inset shows the data at 2 K after 
subtraction of a linear field background and the continuous 
line is the Fourier fit with periods 0.1 T and 0.387 T. These 
periods are independent of temperature within experimental 
resolution. 



If there are fluxons then one should have irreversible 
hysteretic loops of the kind observed in granular super- 
conductors [fj, lD] • Figure [3ja) shows this irreversibility 
that cannot be explained by ferromagnetism, ferroelec- 
tricity due to motion of charges or by usual Abrikosov 
vortices, since no sign of irreversibility has been seen 
within experimental error for magnetic fields applied par- 
allel to the planes. We have a huge anisotropy in an oth- 
erwise a small spin-orbit coupling material. Note that the 
two minima in R are observed at the positive and nega- 
tive fields coming from high fields from the same direc- 
tion. Only by fluxons running between the superconduct- 
ing and the semiconducting-like regions these hysteresis 
loops can be explained. For a better appreciation of the 
hysteresis the inset in Fig. E^a) shows the difference be- 
tween the two curves, i.e. the resistance curve obtained by 
starting at a negative field and sweeping to positive fields 
is subtracted from the resistance curve measured when 
starting at a positive field and sweeping to negative fields. 
The height of the extreme as well as their fields B m (T) 
depend on T. The T-dependence of this irreversibility 
AR as well as B m (T) vanish at Tj ~ 11 K. The reason 
why the irreversible behaviour shown in Fig. [3] vanishes 
at ~ 11 K in contrast to the ~ 25 K observed from the 
oscillatory behavior of Fig. [21 can be easily related to the 
pinning of the fluxons inside the grains. The tempera- 
ture dependence of the irreversibility in field, continuous 
lines in Fig. |3fb), follows (1 — (T/Tj)) 1 ' 5 a similar depen- 
dence as for the irreversibility line of vortices observed in 
high- 1 emp erature superconductors. 

Because it is just graphite, the superconducting re- 
gions have a very small number of free electrons, say 
< 10 -4 electrons per carbon atom [5]. A simple esti- 
mate shows that the London penetration length is larger 
than microns and therefore the Meissner effect should be 
unnoticeable. Also the resistance does not drop to zero 
because the superconducting regions do not percolate, in 
additions to the resistance due to the motion of fluxons. 
The observed hysteresis is a very strong fingerprint of 
superconducting fluxons, difficult to rule out. 

The density of carriers in HOPG samples is very prob- 
ably highly inhomogeneous, and upon region in the sam- 
ple it may be much smaller than 10 _4 /C-atom. What 
might be the physical origin of this superconductivity? 
Graphite contains two carrier families with very differ- 
ent m*, one with a negligible mass called Dirac fermions. 
Therefore, the ratio between masses may be very large, 
100 or larger. These different masses establish large in- 
stabilities if the number of the different carriers is not 
that different. Then we have a large extended band with 
a large Fermi energy corresponding to the light carriers 
and a lower Fermi energy for the heavy carriers, see inset 
in Fig. Q] A large density of the heavy carriers pinned at 
the Fermi energy of the light particles has strong electron 
interaction and creates instabilities that will be discussed 
in other work. In particular for this situation Liu and 
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FIG. 3: (a) Strongly enhanced MR curve near zero field. A 
weak hysteresis appears similar to butterfly MR loops for su- 
perconductors with Josephson-coupled grains Q. For a clear 
observation of the hysteresis we present in the inset the dif- 
ference of the resistance curves (see text), (b) The height of 
the irreversibility maximum AR as well as their field positions 
B m (see inset) vs. temperature. The continuous lines follow 
the function oc (1 - (T/ll)) 1 ' 5 . 



tance up to room temperature, however they are difficult 
to tackle down and their amplitudes change with time, 
an effect that is probably related to the motion of charges 
with current and applied magnetic field. Superconduc- 
tivity in graphite should by no means limited to the 25 K 
here obtained for the small sample, but depends on the 
charge density, defect density and the related instabilities 
at Fermi level. 

We note that hints for superconductivity in HOPG 
samples from SQUID measurements have been invoked 
in the past (T^. However, resolution limits of the mag- 
netometer and the partial admixture of ferromagnetic- 
like signals casted doubts on the origin of those signals. 
Other studies [13| claimed superconductivity in graphite 
based on the metal-insulator transition observed under a 
magnetic field, although superconductivity does not nec- 
essarily need to be invoked to understand this transition. 
There is also a theoretical work that claims high-T c d- 
wave superconductivity in graphite based on resonating 



valence bonds [14 1 



Concluding, in this work we have obtained evidence 
that supports the existence of intrinsic superconductiv- 
ity in HOPG based on the irreversibility of the MR and 
on the quantum oscillations. We think that interior-gap 
- breached superconductivity Q is an interesting start- 
ing concept to understand the observed as well as other 
phenomena in the transport properties of graphite. 
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Humboldt foundation. 



Wilczek |9j have predicted a condensed supcrfluid state 
called interior-gap superconductivity or breached super- 
conductivity. Graphite might be a good candidate where 
some concepts of this theory could be useful. In fact the 
picture they describe for their theory Q is similar to that 
of the inset in Fig. [TJ In this theory no gap exists and 
the material may exhibit p-type superconductivity, which 
has been also discussed for graphite [ll[ as a more robust 
state in a non-homogeneous system. 

The results of Figs. [2] and [3] belong to a micrometer size 
sample (parallel to the planes) and 12 nm thickness in or- 
der to have few potential fluctuations. Measurements in 
two other samples of similar size show similar behavior 
but slightly different T c 's. In larger samples, as for ex- 
ample the other two reported in Fig. 1, the same type 
of effects should be seen but more in terms of universal 
conductance fluctuations. In fact we have observed in 
these and other larger samples fluctuations in the resis- 
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